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Abstract 

We consider the self-similar phase space with reduced fractal dimension d being 
distributed within domain < d < 1 with spectrum f(d). Related thermostatistics 
is shown to be governed by the Tsallis' formalism of the non-extensive statistics, 
where role of the non-additivity parameter plays inverted value f(q) = l/r(q) > 1 of 
the multifractal function r(q) = qd{q) — f(d(q)), being the specific heat, q £ (1, oo) 
is multifractal parameter. In this way, the equipartition law is shown to take place. 
Optimization of the multifractal spectrum f(d) derives the relation between the 
statistical weight and the system complexity. 
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1 Introduction 



Generalization of the thermostatistics is known to be based on the deformation 
procedure of both logarithm and exponential [1,2,3]. The simplest way to 
introduce these functions into the thermostatistics scheme is to consider the 
equation of motion for dimensionless volume 7 = Y /(2nh) 6N of the supported 
phase space (h, N being Dirac-Planck constant and particle number). In the 
course of evolution of the ensemble with statistical weight w = w(j) and 
entropy S, the variation rate of the phase space volume is governed by the 
following equation [3] 



d7 dS 
— = w('y)——. 
dt K u dt 
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From this, one follows the relation dS = Arf/w{ r y) whose integration gives the 
entropy related to the whole statistical weight W: 



j(W) 
7(1) 

Here, we take into account that entropy of a single state S(w = 1) vanishes. 

In the case of the smooth phase space, one has trivial relation w{pj) = 7 
whose insertion into the relation (2) derives to the Boltzmann entropy S = 
InW. However, complex systems has fractal phase space with the dimension 
D < QN, so that relation between the statistical weight and the corresponding 
volume should be replaced by the power-law dependence 

"(7) = iVd (3) 

where reduced fractal dimension d = D/6N < 1 is introduced. Insertion of 
Eq. (3) into integral (2) gives the expression 

S(W)=hi s (dW), WWeI^ (4) 

which is reduced to the Tsallis' logarithm hi q (x) where non-additivity pa- 
rameter q is replaced by inverted value d = 1/d > 1 oi the reduced fractal 
dimension d of the phase space. Naturally, this expression is reduced to the 
Boltzmann entropy in the limit d — > 1. 

Adduced formalism is based on the proposition that the phase space is a 
monofractal determined by single dimension d. The aim of this letter is to 
generalize Tsallis' thermostatistics onto the multifractal phase space with a 
spectrum f(d). Such a generalization for arbitrary distribution f(d) is carried 
out in Section 2. Discussion in Section 3 shows that overall representation of 
the thermostatistics based on the multifractal phase space demands of the 
consistent consideration of input, escort and 'physical' distributions. An opti- 
mization of the spectrum f(d) is considered in Section 4. Section 5 concludes 
our treatment. 



2 Multidractal phase space 

According to the usual recipe [4], the passage to multidractal phase space 
is provided with replacement of the power function (3) by the expression 
ryqd-f(d) whgj-g g [ s multifractal index, multiplier ^"f^ takes into account 
specific number of monofractals with dimension d within the multifractal. As 
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a result, the statistical weight, playing the role of the multifractal measure, 
takes the form 

1 

w q {l) = / l qd - fid) p(d)dd (5) 


where p(d) is the density distribution over dimensions d. Using the method of 
steepest descent, we arrive at the modified power law 

w q ia) ~ Y iq) (6) 

which generalizes the simplest relation (3) due to replacement of the bare 
fractal dimension d by the multifractal function r(q) = qd(q) — f(d(q)). The 
conditions of application of the steepest descent method 



(id 



d 2 f 

= q- 



d=d( q) dd2 



<0 (7) 

d=d(q) 



allow us to find special fractal dimension d(q) related to given parameter q. 

Above consideration shows that the passage from monofractal phase space to 
multifractal one is obtained by replacement of the single dimension d by the 
function r(q) that monotonically increases, taking value r = — 1 at g = and 
t — at q — 1. Limit behaviour of the function r(q) is characterized by the 
asymptotics 

r oc (q - 1) at < q - 1 < 1, r ~ 1 at q — > 00. (8) 

Physical domain of the q parameter variation is bounded by condition q > 1 
which ensures positive values of the function r(q) < 1. 

As a result, we can use well-known Tsallis' formalism of the non-extensive 
statistics where role of the non-additivity parameter plays inverted value 
f(q) = l/r(q) > 1 of the multifractal function r(q) = qd(q) — f(d(q)). Thus, 
the entropy in dependence of the probability distribution Pj has the form [1] 

w i T f-i _ 1 

5 9 = -S^lnT W (Pi), Mx)= __ ■ (9) 
i=i T 1 

With accounting the conditions 

w q w q 

E^ = l, E q = ^ £i P^\ (10) 
i=i i=i 

this expression arrives at the generalized distribution over energy levels 

P i = Z- 1 exp f(q) (-Pe i ) (11) 
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where Z q is the partition function and the deformed exponential is determined 
in the usual manner: 



[1 + (f - l)x]~ at 1 + (f - l)x > 0, 

exp f (x) = <( 

otherwise. 
The normalization condition (10) gives the relation 



(12) 



Wo 



]TPf (<?) = Z-^)- 1 ! - [f(q) - 1] PE q . 



(13) 



i=i 



In spite of the formal similarity of above expressions with the same of the 
Tsallis' statistics, a principle difference takes place: in the last case, the limit 
of extensive systems is reached, when non-additivity parameter takes value 
q = 1, whereas the multifractal phase space is reduced to monofractal one 
with dimension f(q) — > 1 in the opposite limit q — * oo. In this way, the step- 
function spectrum r(q), being r = 1 at q > 1 and r = otherwise, relates to 
the smooth phase space. 

Thermodynamic functions of the model under consideration can be found anal- 
ogously to the Tsallis' non-extensive scheme [1] . However, related expressions 
are very cumbersome even in the simplest case of the ideal gas [5,6,7] and 
amount to the usual form only within the slightly non-extensive limit [8]. At 
the same time, developed scheme allows one to use thermodynamic formalism 
of multifractal objects [9]. Within the latter, the role of a state parameter 
plays the multifractal index q, whose variation may arrive at phase transitions 
if the dependence r(q) has some singularities. It is worthwhile to stress that 
developed scheme arrives directly at related singularities of thermodynamic 
functions type of the internal energy (27) (see below), the entropy (cf. Eq.(4)) 

S q = f(q) \n f(q) (W q ) , f(q) = l/r(q) (14) 

and the free energy 

F q = E q -TS q . (15) 



3 Discussion 



In accordance with the non-extensive thermostatistics [1,5] the distribution 
(11) plays the role of the escort probability related to the input one 



Pi 



p T (<l) _ z~f(q) 



(16) 
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being normalized by the condition 



EpI (9) = 1- (17) 
i=i 



Easily to convince, the expression (16) can be reduced to pseudo-Gibbs form 

Pi = exp r [p g (F g — Ei)} (18) 

if one introduces effective value of the inverted temperature 

(5 q = f(q)Z-^-^ (19) 

and determines the free energy F q in the usual manner: 

F q = -Tln^Z,). (20) 

It is characteristically, the escort distribution (11), determined with the phys- 
ical temperature T = f3~ l , is normalized by the standard condition (10), 
whereas the input probability (16), being deformed exponential (18), is deter- 
mined with non-physical temperature (19) to derive to the deformed free en- 
ergy (20). Thus, both input and escort distributions appear as complementary 
ones. In this way, the escort distribution (11) is deformed with the inverted 
index f(q), whereas the input probability (18) - with the index r{q) itself. 
This exhibits the known 'r(g) <-> f (g) '-duality of Tsallis' thermostatistics [2]. 

According to [5] the physical distribution is neither escort nor input probabil- 
ities, but the following one 



It corresponds to the condition 

w q 

^-E q )P^(e i ) = (22) 
i=i 

instead of the second equation (10). Easily to show, the probabilities (16), 
(21) are connected by the relation 

Vie,) = ZjTM-Mei ~ E q ) (23) 

to derive the distribution 

Vie,) = Z- 1 exp T(g) [-f{q)P ( £i - E q )\ . (24) 

In the case of continuous energy spectrum characterized with the density dis- 
tribution pie), the internal energy related to the condition (22) takes the form 
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E q = J eV{e)p{e)de. (25) 

— oo 

Extreme value of E q is reached at the condition 

P>) V'(e) 



p(e) ~ V(e) 



(26) 



where prime denotes differentiation over e. Usually, the density function is 
reduced to the power law p(e) ~ e cN , c ~ 1, so that p'(e)/p(e) ~ cN/e. Then, 
with using the distribution (24), the condition (26) taken at e = E q arrives at 
the equipartition law 

E q = cr(q)NT (27) 
according to which the value cr(q) is the specific heat. 

Let us stress that above scheme is related to generalized definition of the 
deformed logarithm [2] 

*** )s fw) (28) 

where the function <p{y) is reduced to the statistical weight distribution w q = 
Vi^ fyP e °f Eq.(6). However, this function is obeyed the condition p/ = 
—fiqWqipi) which includes deformed temperature being inconvenient at 
physical considerations. 



4 Optimization of multifractal spectrum 



Up to now, we suppose that the multifractal spectrum f(d) is arbitrary. If it 
is optimized at normalization condition 



/ Mdd 



(29) 



one has to minimize the expression 

7(W q ) 



S q {f(d)} 



7(1) 



1 



qd-f(d) 



p(d)dd 



d7 — 



f(d)dd- 1 



(30) 
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written with accounting Eqs.(2), (5), where E is Lagrange multiplier. As a 
result, we arrive at the equality 

l(W q ) 

J 7 {bd-/(d)]-2r(,) }]ll7d7= ^_ (31) 
7(1) 

whose integration gives, with accounting Eq.(6), the transcendental equation 



\ ^T{q)T q {d)f - [Tq (d) \n(W q ) - 1] - 1 = 0, 

T q (d)=f(q){l-2r(q) + [qd-f(d)]}. 



(32) 
(33) 



This equation is written in the form, when can be used either given spectrum 
function f(d) or the index dependence r(q). In the latter case, we find initially 
the dependence q(d) from the equation 



dr 
dq 



(34) 



q=q(d) 



being conjugated to Eq.(7). Then, substituting this dependence into Eq.(33), 
we arrive at the trivial expression 



T{q) ee T q {d{q)) = f(q) - 1 



(35) 



whose using in Eq.(32) allows us to determine the dependence of the statistical 
weight W q on the complexity E at given function r(q). 

With passage to the smooth phase space, when q — > oo, d — > 1, T(q) — > 0, 
one obtains the statistical weight 



E„o = a^N 



(36) 



which is determined by the specific complexity per one particle. At small 
deviation off the minimum complexity (E — Eoc <^ T,^) and light multifrac- 
tality (1 — r(q) S^, 1 ), linearized equation (32) gives 



W q ~W T{q) \l + r(q) 



1 - 3 (1 - r(q)) E c 



(S-Soo)}, 



w. 



r(q) 



exp<^r(g)E c 



1--(1 

3 V 



r(q)) E c 



(37) 



1 < q < oo. 



In the opposite case r(g) <^ 1, one has with logarithmic accuracy 



W q ~ 



ar(g)E 5 



r(q) 



, a ~ 1, 



(38) 
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5 Conclusion 



As shows above consideration, the thermostatistics of complex systems with 
phase space, whose reduced fractal dimension d is distributed with spectrum 
f(d), is governed by the Tsallis' formalism of the non-extensive statistics. 
In this way, the role of non-additivity parameter plays inverted value of the 
multifractal function r(q) = qd(q) — f(d(q)) which monotonically increases, 
taking value r = at q = 1 and r ~ 1 at q — > oo (the latter limit is related 
to the smooth phase space). The multifractal function r(q) is reduced to the 
specific heat to determine, together with the inverted value f(q) = l/r(q) > 
1, both statistical distributions and thermodynamic functions of the system 
under consideration. At given function r(q), optimization of the normalized 
multifractal spectrum f(d) arrives at the dependence of the statistical weight 
on the system complexity. 
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